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We investigate the production of tp(2S) in nuclear collisions at RHIC and LHC energies. We first 
address charmonium production in 200 GeV d-Au collisions at RHIC; the strong suppression of ip' 
mesons observed in these reactions indicates mechanisms beyond initial cold nuclear matter effects. 
We find that a more complete treatment of hadronic dissociation reactions leads to appreciable ip' 
suppression in the hadronic medium of an expanding fireball background for d-Au collisions. When 
implementing the updated hadronic reaction rates into a fireball for 2.76 TeV Pb-Pb collisions at 
LHC, the regeneration of ip' mesons occurs significantly later than for J/tp’s. Despite a smaller total 
number of regenerated tp ', the stronger radial flow at their time of production induces a marked 
enhancement of their Raa relative to J ftp’s in a transverse-momentum range of p t — 3-6 GeV. We 
explore the consequences and uncertainties of this “sequential regeneration” mechanism on the 
Raa double ratio and find that it can reproduce the trends observed in recent CMS data. 
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I. INTRODUCTION 

Charmonium production in ultra-relativistic heavy-ion 
collisions (URHICs) has been studied for over 30 years. 
The originally proposed J/tp suppression signature of 
Quark-Gluon Plasma (QGP) formation |4] has evolved 
into a more complex problem where both suppression 
and so-called regeneration (or statistical hadronization) 
mechanisms need to be considered. Their interplay and 
relevance depend on collision energy, system size and 
the 3-momentum of the measured charmonia, see, e.g., 
Refs. [HI for recent reviews. The phenomenological 
modeling of these mechanisms, and their relation to the 
underlying in-medium properties, has progressed signif¬ 
icantly in recent years. In particular, kinetic transport 
approaches, when calibrated with existing data from SPS 
and RHIC, have predicted the main features of the J/tp 
production observed in the new energy regime at the 
LHC (although significant uncertainties due to, e.g., 
the open-charm cross section persist @)- These include 
the overall increase of the nuclear modification factor, 
Raa, compared to RHIC energies and its enhancement 
at low transverse momentum, p t . jl] 

Much less is known about the 2 S excited state, 
tp'(36 86). Its small “binding” energy of about 45MeV 
(relative to the DD threshold) renders controlled the¬ 
oretical calculations of its in-medium properties (bind¬ 
ing energy and inelastic reaction rates) challenging. Ex¬ 
perimentally, the tp 1 over J/tp ratio has been measured 
at the SPS 0, where it was found to drop by up 
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to a factor of 3 in central 17.3 GeV Pb-Pb collisions. 
This is consistent with the statistical hadronization ap¬ 
proach Q but it can also be explained by transport 
approaches with lar ge i nelastic reaction rates of the tp’ in 
the hadronic phase Qj), [l3| . More recently, tp ' data have 
become available for 0.2 TeV d-Au collisions at RHIC [l| 
and 5.02 TeV p-Pb collisions at LHC [l|. tp' mesons 
were found to be significantly more suppressed than J/tp 
mesons, which is diffcult to reconcile with initial cold- 
nuclear-matter (CNM) effects since the passing time of 
the highly Lorentz-contracted incoming nuclei is much 
smaller than the formation time scale of the charmonia. 
Consequently, final-state effects have been put forward 
to explain these data, e.g., using the comover suppres¬ 
sion model [0. The latter achieves a good description 
of the collision energy and rapidity dependence of tp' and 
J/tp production in d-Au and p-Pb collisions including ex¬ 
pected shadowing effects on the parton distribution func¬ 
tions (see also Ref. 0 ). 

However, rather unexpected results have emerged from 
recent measurements by the CMS collaboration 0] for 
the double-ratio of the nuclear modification factor, Raa, 
of tp' over J/tp in 2.76 TeV Pb-Pb collisions at the LHC 
(preliminary results are also available from ALICE 0 )- 
At slightly forward rapidities, 1.6<|y|<2.4, and for trans¬ 
verse momenta 3<pt<30 GeV, this double ratio is around 
0.9±0.45±0.3 for semi-central collisions (similar for pe¬ 
ripheral ones), but significantly exceeds one for central 
collisions, 2.3±0.5±0.35. Especially the latter has evaded 
any model explanations thus far, see, e.g., the detailed 
studies in Ref. [20]. On the other hand, around midra¬ 
pidity, and for momenta 6.5<pt<30 GeV, a double ratio 
of around ~0.5 is found, which is much more in line with 
common expectations of a stronger suppression of tp' due 
to its much weaker binding relative to the J/tp. 
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In the present paper we put forward a potential mech¬ 
anism to (partially) resolve the above “puzzle”. Based 
on the rather large inelastic reaction rates for the ip' in 
hadronic matter that we deduce from its suppression in 
d-Au (also in line with the aforementioned SPS data), we 
argue that the inverse reactions of ip' formation in Pb-Pb 
collisions must also happen in the later, hadronic stages 
of the fireball evolution. In particular, the ip' regenera¬ 
tion processes happen later than those for the J/ip whose 
much larger binding energy leads to an earlier “freeze- 
out” than for the ip'. A consequence of such a “sequen¬ 
tial freezeout” is that the collective expansion velocity of 
the medium leads to harder p t spectra for the ip'. Thus, 
in terms of the Raa , the ip' can outshine the J/ip in a 
momentum range of pt > M,/,, which happens to coincide 
with the lower CMS pt cut. On the other hand, at higher 
p t , the regeneration contribution ceases giving way to a 
more sequential-like suppression pattern of primordially 
produced charmonia. 

Our paper is organized as follow: In Sec. El we re¬ 
visit our hadronic reaction rates for both ip' and J/ip, 
apply them within a schematic fireball for d-Au collsions 
at RHIC and compare our results to PHENIX data. In 
Sec. ED we implement the updated hadronic rates in our 
earlier constructed transport approach. In particular, we 
estimate the time windows for sequential charmonium re¬ 
generation, and elaborate the uncertainty for the under¬ 
lying charmonium p t spectra in the context of the CMS 
data for the ip'/ J/ip Raa double ratio. We conclude in 
Sec. m 


II. HADRONIC DISSOCIATION OF 
CHARMONIA 

Dissociation rates of charmonia in hadronic matter are 
usually considered to be much smaller than in the QGP 
(see, e.g., the discussion in Ref. i)- However, for the 
ip' this is not so obvious, since the proximity of its mass 
to the DD threshold provides a large phase space for 
break-up reactions. In the following, we revisit hadronic 
reactions rates for J/ip and ip' mesons based on effective 
meson Lagrangians (Sec. Ill AD and evaluate their conse¬ 
quences for final-state effects in d-Au reactions at RHIC 

(Sec. ESI)- 

A. Update on Dissociation Rates 

Our startingpoint is the previously employed hadronic 
reaction rate [131 based on flavor-S't/(4) meson exchange 
models [El [ 22 J for the processes J/ip + p —> D + D 1 D* + 
D* (exothermic for mj+ m p > mn + mp, and en¬ 
dothermic for mj/^ + m p < Tod* + mp,*) and J/ip + n -A 
D* + D, D + D* (endothermic). For those reactions, 
the J/ip dissociation rate at T=170MeV amounts to 1- 
2 MeV, corresponding to a lifetime of 100-200 fm/c. Even 
with the typical uncertainties (factor ~2-3) associated 
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FIG. 1: (Color online) Temperature dependence of hadronic 
dissociation rates for J/ip (upper panel) and ip' (lower panel) 
at rest in a thermal bath. Previous results for a up gas with 
A=lGeV (solid red lines) are compared to our updated re¬ 
sults for a meson resonance gas using A=1 GeV and 2 GeV 
(solid green and blue lines, respectively). The blue-dotted 
lines additionally account for finite meson chemical potential 
that build up for temperatures below the chemical freezeout 
in URHICs 0. 


with the hadronic formfactor cutoff values, this is too 
small to affect the J/ip abundance during the 5-10 fm/c 
lifetime of the hadronic phase in URHICs. For the ip ', ge¬ 
ometric scaling by the vacuum radius has been assumed, 
increasing its rates by a factor (r^/ /rj^,) 2 , which is ap¬ 
proximately compatible with constitutent-quark model 
calculations [El . 

A hadron resonance gas (HRG), however, contains 
many more species than 7 r and p. To estimate their 
impact on the charmonium dissociation rates, we sim¬ 
ply adopt the existing p- and 7r-induced matrix elements 
and shift their kinematics according to the pertinent 2- 
particle threshold, i.e., 

rfh/4 T ) = J (0/- Y (^(fc);T)< +J/ 4 s , St Y hr )u rel 

(i) 
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FIG. 2: (Color online) Temperature evolution in 0.2 TeV d- 
Au collisions assuming a thermal fireball model. The blue 
(red) curve is for an impact parameter of 6=0(6) fm. 


with s = {pj/^j + k) 2 and s^ r = (m J/v , + m x ) 2 for 
exothermic and s^ r = (2mu) 2 for endothermic channels 
(s(b r = ( mo + w_d s ) 2 if X contains a strange quark). 
We include a total of 52 non-strange and single-strange 
meson species, up to a mass of mx = 2GeV. As be¬ 
fore, we apply geometric scaling to obtain the reaction 
rates for the \c and ip'. Our results for a meson gas in 
chemical equilibrium are summarized by the solid lines 
in Fig.[T| At the highest temperature (T=180MeV), the 
additional resonances enhance the J/ip dissociation rate 
by a factor of ~2.5, and another factor of ~2.5 when in¬ 
creasing the hadronic formfactor cutoff from A=1 GeV to 
2 GeV, reaching a maximal rate of 10.5 MeV. For the ip', 
geometric scaling leads a maximun rate of up to 35 MeV, 
translating into a lifetime of ~6fm/c which is now com¬ 
parable to the duration of the hadronic phase in URHICs. 
This becomes even more significant if chemcial freezeout 
is accounted for, which implies the build-up of meson 
chemical potentials leading to a slower decrease of the 
meson densities as temperature decreases, cf. dotted lines 
in Fig. [□ 


B. Charmonium Production in 0.2 TeV d-Au 
Collisions 


We proceed by implementing our updated hadronic re¬ 
action rates into the thermal-rate equation framework de¬ 
veloped in Ref. 25J. Another new aspect relative to our 
previous work [26] is that we allow for final-state effects 
in small collision systems, specifically for d-Au at RHIC, 
where a rather large ip' suppression has been observed 
(while J/ip’s are much less suppressed). 

Toward this end, we construct a thermal fireball us¬ 
ing the same methods as before for AA collisions. We 
determine the total entropy of the fireball by matching 
the observed final-state hadron abundancies. Based on 


FIG. 3: (Color online) Centrality dependence of the nuclear 
modification factor for charmonia in 0.2 TeV d-Au collisions 
at RHIC; CNM effects for all charmonia are shown by the 
black dotted line (consisent with Ref. [Til ] based on the EPS09 
shadowing [Til] parameterization); the blue band illustrates 
uncertainties in both QGP and hadronic dissociation rates of 
the ip ', while the red band illustrates the uncertainties due 
to the hadronic dissociation of the J/ip (solid (dotted) line: 
A=0.75(2) GeV). The PHENIX data are from Refs. [H|,[3(|. 


a Glauber model for the centrality dependent nuclear 
overlap function, we initialize the fireball with a trans¬ 
verse radius of i?o—2.5fm/c. The expansion is then mod¬ 
eled with the same acceleration as in 0.2 TeV Au-Au be¬ 
fore [26}. For simplicity, we stick with our first-order 
equation of state, transitioning from a quasi-particle 
QGP into a HRG through a mixed phase at T= 180 MeV 
(we do not expect large changes when utilizing a modern 
cross-over EoS, which has been checked for dilepton [27J 
and bottomonium observables [Hj). Kinetic freezeout is 
also constructed as before, with a freezeout temperature 
mildly decreasing with centrality, from Tf 0 ~155MeV in 
peripheral to ~142MeV in central collisions, resulting in 
fireball lifetimes of 0.5-3 fm/c, cf. Fig. [2] For the most 
central collisions, a short QGP phase with initial temper¬ 
ature To=190MeV is followed by a lfm/c mixed phase 
and a 1.5 fm/c hadronic phase. 

The initial CNM effects are assumed to be identical 
for all charmonia, essentially given by a shadowing sup¬ 
pression fl6l ] which we mimic with an “effective” nuclear 
absorption cross section of <7 a b s =2.4mb. Our results for 
the centrality dependence of the nuclear modification fac¬ 
tor, 


^AA^part) 


Nqf(N pa rt) 

N coll (-N part) ’ 


( 2 ) 


are summarized for both T = J/ip,ip' in Fig. [Sj in com¬ 
parison to PHENIX data [h] (IV co n: number of primor¬ 
dial NN collisions). For the J/ip, there is a moderate 
suppression beyond CNM effects; the additional final- 
state effects are due to a small QGP suppression on the 
direct J/ip's, as well as a suppression of the feeddown 
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from Xc’s and ip', with little room for additional hadronic 
suppression (thus, in principle, preferring a small value 
for the hadronic cutoff parameter). On the other hand, 
for the ip', our baseline QGP+HRG suppression is not 
enough to account for the marked suppression beyond 
CNM effects. Here, a large formfactor cutoff of A=2 GeV 
is preferred to augment the hadronic suppression. An 
additional increase of the QGP suppression rate of the 
ip' by a factor of 2 could also be helpful (such an in¬ 
crease may arise, e.g., from nonperturbative heavy-quark 
interactions with light partons). Further increasing the 
hadronic rate by a factor of 2 1 has a similar effect, re¬ 
ducing the discrepancy with the most central PHENIX 
datum. 

We finally note that a hadronic ip' dissociation rate 
well beyond the one from the np gas (by a factor of 5 
or more) was previously deduced within our setup [l3j 
to be able to account for the ip' suppression observed in 
S-U and Pb-Pb collisions at the SPS [10]. Our newly 
calculated rates in the present ms. are in line with this 
notion. 


III. SEQUENTIAL REGENERATION OF 
CHARMONIA 

In this section, we investigate the consequences of the 
updated hadronic reaction rates within our previously 
constructed thermal fireball expansion (Sec. IIII All , fol¬ 
lowed by a more generic evaluation of the associated un¬ 
certainties specifically in the context of the ip' / J/ip Raa 
double ratio ('Sec. lIII 51 ) . We focus on 2.76 TeV Pb-Pb col¬ 
lisions at the LHC. Our earlier predictions for inclusive 
J/ip production in these reactions [5] resulted in a fair 
agreement with the centrality, transverse-momentum and 
rapidiy dependencies observed by the ALICE and CMS 
collaborations, and thus serves as our framework to eval¬ 
uate ip' observables. Since the pertinent CMS data are 
for “prompt” J/ip and ip' production, we do not include 
contributions from B feeddown. 


A. Fireball Model 

For our fireball results we focus on the so-called 
“strong-binding scenario”, where the in-medium charrno- 
nium properties are taken with guidance from a T-matrix 
approach [3l| with the internal energy from lattice-QCD 
as underlying potential. This assumption gives a better 
agreement than using the free energy both with correla¬ 
tors from lattice-QCD and the overall charmonium phe¬ 
nomenology at SPS and RHIC {Hj. For definiteness, we 


1 This could be due to inelastic reactions with baryons and an- 
tibaryons or direct ij>‘ —> DD decays with an in-medium reduced 
D meson mass [2. r ( . neither of which we have calculated here. 



time(fm/c) 

FIG. 4: (Color online) Time dependence of J/ip and i />’ nu¬ 
clear modification factors in central(A pa rt=324) Pb-Pb col¬ 
lisions at -ys=2.76TeV (the denominator for the J/ip Raa 
includes 40% feeddown from excited states while the numer¬ 
ator does not). The red (blue) curves are for direct J/ip {ip'), 
with the solid (dotted) line styles representing the regener¬ 
ation (primordial) contributions; the dash-dotted curves in¬ 
dicate the pertinent equilibrium limits (including a thermal 
relaxation time correction 0), starting from the time when 
the fireball has cooled down to the dissociation temperature 
below which regeneration commences. 


employ the hadronic rates with A=0.75(2) GeV for the 
J/ip {ip'), and a factor of 2 increased QGP rate of the ip'. 

Let us first inspect the time evolution of direct J/ip 
and ip' mesons in 0-20% central Pb-Pb at midrapidity 
(without shadowing), as following from the solution of 
the kinetic rate equation, 

^ = -r*[IV*(T) - iV; q (T(r))] , (3) 

see Fig. [4] Compared to our previous results (cf. lower 
panel in Fig. 2 of Ref. 1), the J/ip now picks up a re¬ 
generation contribution in the hadronic phase, by about 
0.15 units in Raa, which, despite a small rate, is due 
to the large equilibrium limit. Most of the production, 
however, occurs prior to the onset of the mixed phase 
at r ~ 5.5fm/c. On the contrary, ip' production only 
starts to set in at that point, leveling off at around 9- 
10 fm/c, when the temperature of the fireball has dropped 
to about 150-160 MeV. The main qualitative and robust 
feature here is that the lower dissociation temperature 
of the ip' , relative to the J/ip, implies a later production 
through regeneration in the time evolution of the fireball 
in URHICs. 

The sequential regeneration of J/ip and ip' has rather 
dramatic consequences on their transverse-momentum 
{pt) spectra. Following our previous work ffj, we ap¬ 
proximate the pt spectra of the regeneration components 
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with the standard blast-wave expression, 

Pt sinh p(r) 
T ' 

(4) 

implying thermalized charm-quark distributions. The 
rrit = y/pt + m 2 is the transverse mass. As our default, 
we evaluate this expression when most of the pertinent 
charmonium yield has built up, i.e., at T c for the J/ip 
(r=5.5fm/c in Fig. [4j and at T=155MeV for the ip' 
(r=9.7fm/c in Fig. (U). The resulting nuclear modifi¬ 
cation factors, Raa( pt ) , for both regeneration and sur¬ 
viving primordial components are displayed in Fig. [5] 
(we have assumed the same initial spectra for J/ij) and 
ip' from pp collisions, figuring into the denominator of 
Raa)- Clearly, there is a significant uncertainty associ¬ 
ated with this procedure which we address in Sec. IIII Bl 
below (we recall, however, that our pertinent predictions 
for the J/ip gave fair agreement with the observed pt 
spectra at LHC). Three main qualitative features can be 
gleaned from comparing the Raa’s for J/ip and ip' in 
Fig. [5] (a) At low < 3 GeV, the regeneration yield 
of the J/ip dominates over the one from ip 1 , as a conse¬ 
quence of the approach toward equilibrium which favors 
the smaller J/ip mass; (b) At intermediate pt — 3—6 GeV, 
the significantly harder blast wave for the ip' generates 
a shift of the “flow bump” which exceeds the regener¬ 
ation contribution in the J/ip Raa', (c) at still higher 
momenta, > 6 — 8 GeV, regeneration gives way to 
(suppressed) primordial production, where the J/ip Raa 
exceeds again the one from the more weakly bound ip'. 
Items (b) and (c) are in qualitative agreement with the 
trends observed in the pertinent Raa/^-aa double ratio 
observed by CMS [3. In the next section we explore 
some of the uncertainties in our calculations. 


B. Schematic Model 

To better quantify variations in the interplay of the dif¬ 
ferent production components of both J/ip and ip' , let us 
first formulate a baseline scenario motivated by existing 
experimental data for the J/ip and the thermal fireball 
calculations in the previous sections. For simplicity, we 
assume the primordial parts to be constant in R\A(pt)- 
This reflects our currently limited knowledge about the pt 
dependence of the dissociation rates (see, e.g., Ref. 0), 
formation time effects, etc. For 0-20% Pb-Pb collision we 
take 0.15-0.25 for the J/ip (compatible with high-p t CMS 
data for prompt J/ip [32|) and 0-0.075 for the ip ', to re¬ 
flect its stronger absoprtion as a loosely bound state. For 
the regeneration components, we choose total yields such 
that the total momentum-integrated Raa amounts to 
0.55-0.65 for the J/ip (compatible with ALICE data (HI) 
and 0.4 for the ip' (as suggested by our fireball results). 
For the sequential freezeout, which determines the tem¬ 
perature and flow strength in the blast-wave spectra of 
the regeneration components, we employ the correlation 



Pt(GeV) 


FIG. 5: (Color online) Nuclear modification factor for J/ip 
(red lines) and ip' mesons (blue lines) as a function of pt for 0- 
20% Pb-Pb(2.76 TeV) collisions. The total Raa for each me¬ 
son (solid lines) is decomposed in a regeneration component 
(dotted lines, evaluated with a blast-wave ansatz in the fire¬ 
ball evolution at sequential freezeout times) and a suppressed 
primordial component (dash-dotted lines, as obtained from a 
Boltzmann transport equation without gain term), cf. Ref. [s| 
for further details. 



Pt(GeV) 


FIG. 6: (Color online) Transverse-momentum spectra of J/ip 
(red lines) and ip' mesons (blue lines) for both suppressed 
primordial (dashed lines, assuming Raa =0-2 and 7 ?a A = 0.05) 
and regeneration components (solid lines); the total yields 
(not the shapes) are normalized to their A co n-scaled number 
in pp collisions. 


given by the expanding fireball, i.e., r ~ 5.5(9.7) fm/c for 
the J/ip (ip') corresponding to Tf o =180(155) MeV and 
F=0.32(0.46)c. Recall that the effective slope parame¬ 
ter for the blast-wave spectra is approximately given by 
T e ff « T + m^v 2 , which is mostly driven by the flow term 
due to the large mass of the charmonia. In Fig. [G] we dis¬ 
play the pertinent pt spectra for central Pb-Pb(2.76 TeV), 
normalized by the V co n-scaled yields in pp , to mimic the 
relative magnitude of ip' and J/ip contributions in their 


rfivr 

Ptdpt 


= No(b)mt / rdrKi( 


m t cosh p(r) 
T 


)M 
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FIG. 7: (Color online) Ratio of the Raa for ip 1 to the one for 
J/ip (the so-called “double ratio”) as a function of centrality 
in 2.76 TeV Pb-Pb collisions. The red and blue boxes are our 
schematic-model results for the double ratio in the momen¬ 
tum range p t =3-30 GeV and pt= 6.5-30 GeV, respectively, and 
are compared to CMS data [l8|. In the upper (lower) panel, 
the boxes indicate theoretical uncertainties in the modeling of 
the primordial (regeneration) component of the charmonium 
yields and spectra. The vertical purple lines represent our 
“best fit” in the upper panel, and a “realistic” baseline for 
the blast-wave variations in the lower panel. 


respective Raa’s . The plot highlights again the main ef¬ 
fect proposed in this paper: due to the stronger flow for 
the regenerated ip', its Raa can rise above the one for 
the J/ip in a limited pt window around 3-6 GeV, before 
primordial production takes over again. 

For semi-central Pb-Pb (20-40%) we construct the 
baseline following the same reasoning as for central colli¬ 
sions. We increase the primordial Raa’s for J/ip to 0.35- 
0.45 [Hj and to 0.1-0.2 for ip' , and evaluate the blast- 
wave pt spectra at the same freezeout temperatures as 
for central collisions (but with flow velocities given by 
the fireball expansion for 20-40% centrality). 

The resulting double ratios covering the above- 
specified ranges in the primordial components are shown 
in the upper panel of Fig. [7] and compared to CMS 


data [HI (we did not include the ALICE data [l9[ as 
they contain feeddown contributions from 13-meson de¬ 
cays, which could become rather significant at high p t , 
see, e.g., Ref. [2()j for a calculation including those). The 
basic trends of the data can be reproduced within our ap¬ 
proach. For defintieness, we quote (approximate) “best 
fit” values of 0.15 and ~0 (0.35 and 0.2) for the primor¬ 
dial J/ip and ip’’ values, respectively, in (semi-) central 
Pb-Pb collisions, resulting in the purple horizontal bars 
in the upper panel of Fig. [71 

Finally, we illustrate the sensitivity of the double ratios 
to the blast-wave parameters, by varying the freezeout 
temperatures over the ranges T f ^^=180-200MeV and 
Tf o =150-165 MeV, along with the pertinent flow veloc¬ 
ities from the fireball model, and fixing the primordial 
components at Raa =0.2(0.35) and R^a =0.05(0.15) for 
(semi-) central Pb-Pb. The larger regeneration compo¬ 
nents and larger flow velocities in central collisions render 
the double ratios more sensitive to the details of the se¬ 
quential regeneration mechanism. 


IV. CONCLUSIONS 

In the present work, we have investigated the produc¬ 
tion systematics of ip' mesons in URHICs. We first revis¬ 
ited the problem of hadronic ip' dissociation and found 
that a more complete inclusion of hadronic states in a 
resonance gas suggests a marked increase of its inelastic 
reaction rates. When implementing these rates into an 
expanding fireball for d-Au collisions at RHIC, we found 
a much improved description of the rather strong sup¬ 
pression of ip' mesons observed in these reactions. This 
is similar in spirit to, and thus supports, the recently 
suggested comover suppression effects [3 i n dA and pA 
reactions at RHIC and LHC. We then evaluated ip' trans¬ 
port in Pb-Pb collisions at the LHC using our existing 
fireball approach which previously provided fair agree¬ 
ment with J/ip data. The key features in our approach 
are the lower dissociation temperature of the ip' relative 
to the J/ip and its sizable hadronic reaction rates. This 
implies a sequential freezeout of those two mesons, with 
most of the ip' regeneration occurring later in the fire¬ 
ball evolution. The larger collective medium flow then 
leads to an enhancement of the ip' regeneration yield in a 
Pt region around 3-6 GeV, transitioning to (suppressed) 
primordial production at higher p t . While quantitative 
predictions of this mechanism are beyond current the¬ 
oretical control, we have shown that variations in the 
ingredients to the sequential regeneration scenario pro¬ 
duce trends in the ip' -over- J/ip Raa double ratio which 
agree with recent CMS data. We therefore believe that 
the qualitative features of this mechanism are robust and 
provide a candidate to contribute to the understanding of 
these data. In fact, if corroborated, sequential regenera¬ 
tion may serve as a tool to extract in-medium properties 
of the ip' from URHIC data. 
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